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Abstract 

Some introduced populations thrive and evolve despite the presumed loss of 
diversity at introduction. We aimed to quantify the amount of genetic diversity 
retained at introduction in species that have shown evidence of adaptation to 
their introduced environments. Samples were taken from native and introduced 
ranges of Arctotheca populifolia and Petrorhagia nanteuilii. Using microsatellite 
data, we identified the source for each introduction, estimated genetic diversity 
in native and introduced populations, and calculated the amount of diversity 
retained in introduced populations. These values were compared to those from 
a literature review of diversity in native, confamilial populations and to esti- 
mates of genetic diversity retained at introduction. Gene diversity in the native 
range of both species was significantly lower than for confamilials. We found 
that, on average, introduced populations showing evidence of adaptation to 
their new environments retained 81% of the genetic diversity from the native 
range. Introduced populations of P. nanteuilii had higher genetic diversity than 
found in the native source populations, whereas introduced populations of 
A. populifolia retained only 14% of its native diversity in one introduction and 
1% in another. Our literature review has shown that most introductions dem- 
onstrating adaptive ability have lost diversity upon introduction. The two spe- 
cies studied here had exceptionally low native range genetic diversity. Further, 
the two introductions of A. populifolia represent the largest percentage loss of 
genetic diversity in a species showing evidence of substantial morphological 
change in the introduced range. While high genetic diversity may increase the 
likelihood of invasion success, the species examined here adapted to their new 
environments with very little neutral genetic diversity. This finding suggests that 
even introductions founded by small numbers of individuals have the potential 
to become invasive. 



Introduction 

Biological invasions present one of the greatest environ- 
mental challenges of our time, yet the drivers of successful 
invasion remain poorly understood. The concept that 
genetic diversity in the founding population is positively 
related to the probability of invasion success is one of the 
oldest hypotheses in invasion biology (e.g., Mayr 1965) 
and continues to be supported by recent research 



(Crawford and Whitney 2010; Jones and Gomulkiewicz 
2012). However, the debate over the importance of genetic 
diversity to invasion success continues (Sakai et al. 2001; 
Kolbe et al. 2004; Roman and Darling 2007; Hufbauer 
2008). This hypothesis presents a paradox because genetic 
bottlenecks are expected to occur at introduction, reduc- 
ing the potential for introduced populations to adapt to 
novel environments (Allendorf and Lundquist 2003; Allen- 
dorf and Luikart 2007), but despite this, many introduced 
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populations thrive. In some cases, this has been explained 
by high levels of propagule pressure through multiple 
introduction events, resulting in introduced populations 
having greater genetic diversity than is found in the native 
range (Kolbe et al. 2004; Genton et al. 2005). However, 
successful invasions are not always accompanied by high 
genetic diversity and sometimes are depauperate in neutral 
genetic variation (Ren et al. 2005; Mergeay et al. 2006; 
Zimmermann et al. 2010). 

Many plant and animal populations expanding into 
novel environments not only thrive but also exhibit rapid 
evolutionary changes in crucial traits such as dispersal 
ability, reproductive output, phenotypic plasticity, and 
size (Blossey and Notzold 1995; Cody and Overton 1996; 
Siemann and Rogers 2001; Bossdorf et al. 2005; Phillips 
et al. 2006; Richards et al. 2006; Cheptou et al. 2008; Rid- 
ley and Ellstrand 2009; Buswell et al. 2011). This empiri- 
cal evidence is supported by simulations, demonstrating 
that evolution may move at a faster rate when an organ- 
ism's environment varies (Kashtan et al. 2007) and inva- 
sive populations often experience extreme environmental 
shifts. Further, it appears that rapid evolution in invasive 
species may be quite common. For example, Buswell et al. 
(2011) studied herbarium specimens of 23 plant species 
introduced to Australia and sampled repeatedly across the 
past -150 years to identify evidence of significant 
morphological change. They concluded that changes had 
occurred in 70% of these species following their introduc- 
tion and that this was most likely the result of rapid 
evolution. Evidence of rapid evolution in novel environ- 
ments supports the idea that genetic diversity is impor- 
tant to the success of introduced populations because 
adaptations following introduction are more likely to be 
derived from standing genetic variation rather than muta- 
tion (Barrett and Schluter 2008). Nevertheless, several 
studies have demonstrated rapid evolution in the presence 
of low genetic diversity in introduced ranges (Dlugosch 
and Parker 2008b; Harris et al. 2012), suggesting that the 
level of standing genetic diversity required for adaptation 
may, in fact, be quite low. 

Quantitative genetic theory predicts that the extent of 
adaptive genetic change due to pre-existing genetic varia- 
tion in the initial population in a new environment, as 
well as adaptation due to new mutations arising in the 
new environment, will be an increasing function of selec- 
tion, genetic diversity, genetically effective population size, 
and number of generations (Robertson 1960; Weber 2004). 
The extent of adaptation expected for both pre-existing 
diversity and novel mutation can be predicted (see Appen- 
dix SI for details), and these predictions are supported by 
empirical evidence (Frankham 1980b, 1983; Weber 2004; 
Frankham et al. 2010). As these predictions assume that 
genetic variation is neutral, genetic adaptation should 



increase with levels of neutral genetic variation, other fac- 
tors being equal (Frankham et al. 1999, 2010). 

In this study, we examine the relationship between 
neutral genetic diversity and rapid evolution in intro- 
duced species, using two species that have exhibited 
significant morphological change since their introduction 
to Australia. First, unlike many other studies of rapid 
evolution after introduction, we aimed to determine the 
exact source population(s) for the introductions for both 
species. This information is important for the accurate 
comparison of the genetic characteristics of the intro- 
duced populations from the actual source populations in 
their native ranges. Second, we characterized genetic 
diversity in the native and introduced ranges of both 
species. Finally, we surveyed the literature to determine 
( 1 ) whether the levels of genetic diversity we found in the 
native range of both species were typical of other species 
within those families and (2) whether our species retained 
similar levels of genetic diversity at introduction com- 
pared with other introduced species showing evidence of 
rapid evolution in their introduced environments. We 
expected that genetic diversity in the species studied here 
would not be low relative to their families because they 
had demonstrated the ability to undergo morphological 
change in their introduced environments. We also 
expected the change in diversity between the native and 
introduced populations in our study species to be similar 
to that of other introduced species showing evidence of 
rapid evolution. 

Methods 

Study species 

We aimed to select species that had shown potential for 
rapid evolution through postintroduction morphological 
change. We chose species with restricted native and intro- 
duced ranges so that we could comprehensively sample 
across their distributions. Annual or short-lived perennial 
species with sexual reproduction were selected, because 
these species have had more generations since introduc- 
tion, increasing the opportunity for evolution to occur in 
the introduced range. We avoided selecting crop and pas- 
ture species that were likely to have been introduced 
many times. 

Arctotheca populifolia (Fig. 1), chosen based on the 
findings of rapid morphological change in introduced 
populations (Buswell et al. 2011), is in the Asteraceae and 
is a perennial, herbaceous succulent native to South 
Africa and introduced to Australia (Harden 1992). The 
first records of this species in Australia occurred in the 
1930s on both the east and west coasts (AVH Database, 
2012). The Australian distribution encompasses coastal 
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Figure 1. Australian sample of Arctotheca populifolia (photograph by C. Brandenburger). 



environments from Geraldton (Western Australia) to 
northern New South Wales (Heyligers 1998). The Global 
Compendium of Weeds lists A. populifolia as an agricul- 
tural and environmental weed (GCW database 2012). 

The second study species was Petrorhagia nanteuilii, 
which is an annual, herbaceous plant in the Caryophylla- 
ceae. It is native to western Europe and western North 
Africa (Ball and Heywood 1964) and introduced to Asia, 
Australia, North America, South America, and Macarone- 
sia (GRIN database 2012). This species was first recorded 
in Australia in 1882, and, currently, the Australian distri- 
bution is restricted to the southeast, ranging from Bris- 
bane to Adelaide (AVH Database, 2012). Petrorhagia 
nanteuilii is also listed in the Global Compendium of 
Weeds as an agricultural and environmental weed (GCW 
database 2012). 



Similar to A. populifolia, P. nanteuilii also showed 
evidence of morphological change over time since intro- 
duction (see Results, below). This was determined using 
herbarium specimens following methods described by Bu- 
swell et al. (2011). We measured height on all available 
specimens at the National Herbarium of Victoria (MEL) 
at the Royal Botanic Gardens, Melbourne. This gave data 
for 184 plants from 56 herbarium sheets ranging in collec- 
tion date from 1882 to 1998. No leaf traits were measured 
because leaves do not press well in this species. All plants 
had grown in the range of the introduction, in Victoria 
and New South Wales in Australia. We ran a general linear 
model including region and year as predictors and log 10 - 
transformed height as a dependent variable. The term for 
region was included to prevent the possibility that a popu- 
lation expansion along an environmental gradient would 
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Figure 2. Sampled areas with place name 
abbreviations and number of individuals 
sampled in parentheses. (A) native range 
samples of Arctotheca populifolia. (B) 
Australian introduced range samples of 
A. populifolia (triangles) and Petrorhagia 
nanteuilii (squares). (C) native range samples of 
P. nanteuilii. Genetic groups are indicated by 
bars labeled with group name (i.e., A,; see 
Results). Note that group assignment of P. 
nanteuilii samples PM and Gl is ambiguous 
(see Figs 3 and 4). 



be mistaken for adaption to the native range across time 
(Buswell et al. 2011). To do this, we recorded the region of 
origin for each sample. Because most regions were repre- 
sented by relatively few specimens, we pooled bioregions 
to construct four broad climate regions: (1) humid coast 
and hinterlands (including East Gippsland, Victoria, and 
the New South Wales Central Coast and South Coast), (2) 
humid highlands (including Eastern Highlands, the Snow- 
fields, and the Southern Tablelands), (3) subhumid slopes 
(including the Victorian Midlands and Riverina, and the 
New South Wales South West Plains, South West Slopes, 
and North West Slopes), and (4) semi-Mediterranean 
(including the Victorian Volcanic Plain, the Grampians 



and Wannon). In order to acknowledge the nonindepen- 
dence of plants from the same herbarium sheet, we 
weighted individuals according to the number measured 
on the herbarium sheet such that the weights for all the 
plants on each sheet sum to one. For example, a single 
plant on a sheet received a weight of one, while two indi- 
viduals on the same sheet each received a weight of 0.5. 
Analyses were performed in JMP, version 5 (SAS Institute, 
Cary, NC). 

For both species, we also measured plant height over 
time in the native range, using the methods described 
above. This was done in order to determine whether 
any changes identified in the introduced range were 
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Figure 3. Log 10 plant height of Petrorhagia nanteuilii introduced to Australia measured from herbarium specimens sampled from 1880 to 2000, 
classed by climatic region. Values increased significantly across time (weighted general linear model including a term for region; R 2 = 0.06; 
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concurrently occurring in the native range, perhaps as a 
result of global climate change. For these data, region was 
not included as a term due to the small number of sam- 
ples available for each region. In total, 52 herbarium sam- 
ples from 28 sheets were measured from the native range 
of A. populifolia and 86 samples from 26 sheets for native 
range P. nanteuilii. 



Genetic sampling 

We sampled leaves from 348 A. populifolia plants from 10 
sites covering the native range (N = 188; Fig. 2A) and 
seven sites across the introduced range in Australia 
(N = 160; Fig. 2B, triangles). For P. nanteuilii, we sam- 
pled a total of 345 plants, including those from 12 sites 
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Figure 4. Structure analyses. Evanno et al.'s (2005) AK values for each putative number of populations (K) for (A) Arctotheca populifolia and (B) 
Petrorhagia nanteuilii. Structure Q plots generated using the maximum value of AK indicate A. populifolia (C) and P. nanteuilii (D) samples 
represent two genetic groups each, demarcated by labeled bars (i.e., A,). Samples from the introduced range are denoted by asterisks. Each 
individual is represented by a vertical line showing degree of admixture. Sample name abbreviations are defined in Figure 1. 
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in the native range (N = 282; Fig. 2C) and two sites 
across the introduced range in Australia (Fig. 2B, 
squares). Two attempts were made to sample this species 
in the vicinity of Adelaide, South Australia, at the west- 
ernmost reported extreme of the Australian distribution, 
and in the vicinity of Sydney where P. nanteuilii has also 
been reported; however, on all occasions, none were pres- 
ent. Leaves were placed in vials containing a solution of 
40% sodium chloride, 4% sodium ascorbate, 4% silica, 
and 3% cetyltrimethylammonium bromide (Thompson 
2002) and stored at 4°C. To prepare samples for extrac- 
tion, leaves were removed from the preservative, washed 
in Milli-Q water, patted dry, and frozen at — 70°C prior 
to freeze drying. Freeze-dried samples were crushed and 
DNA was extracted using a NucleoSpin 96 Extraction II 
Kit (Macherey-Nagel, Diiren, Germany). 

Microsatellites were developed using next-generation 
sequencing on the GS-FLX 454 platform (Roche, Manheim, 
Germany) following methods described by Abdelkrim et al. 
(2009). QDD v 0.9.0.0 Beta (Meglecz et al. 2010) was used 
to identify microsatellites, and primers were designed using 
the program Primer 3 (Rozen and Skaletsky 2000). A panel 
of polymorphic markers was chosen for each species (A. 
populifolia, seven microsatellite loci; P. nanteuilii, 12 micro- 
satellite loci; Table SI). Using universal primers (Neilan 
et al. 1997) having four differently colored fluorescent 
labels, we multiplexed PCRs within label color and 
multiloaded all loci for each species into a single reaction 
per individual. The step-down PCR protocol consisted of 
ten cycles each at the following annealing temperatures: 
70°C, 64°C, 58°C, 54°C, 50°C. Samples were genotyped 
using an ABI 3730 (Applied Biosystems, Foster City, CA) 
using GS-500 (Liz) in each capillary as a size standard. 
Allele sizes were estimated on GeneMapper, version 3.7 
(Applied Biosystems). 

Statistical analyses of genetic data 

We tested microsatellite data for departures from Hardy- 
Weinberg and linkage equilibrium in Arlequin, version 
3.5.1.2 (Excoffier et al. 2005), and P-values were Bonfer- 
roni corrected. We used Structure, version 2.2 (Pritchard 
et al. 2000; Falush et al. 2003), to determine whether mul- 
tiple genetic groups were present across the range of each 
species and to determine the native source of introduced 
populations. For this analysis, we used the admixture 
model with correlated allele frequencies and tested the 
number of genetic groups (K) for each value of K between 
one and ten. We ran ten replicates for each value of K, 
each run having a burn-in period of 100,000 Markov chain 
Monte Carlo steps followed by 10 6 iterations. The most 
likely number of genetic groups was inferred using Evanno 
et al.'s (2005) AK method. We determined group member- 



ship assignment of each sample using the highest propor- 
tion of membership across all ten runs of Structure. 
Principal coordinate analysis (PCoA) conducted in GenA- 
lEx v. 6.3 (Peakall and Smouse 2006) was used to visualize 
genetic distances (Nei 1972) between populations. 

Many authors have stressed that a spectrum of diversity 
measures gives the best summary of diversity (Pielou 
1966; Hill 1973). Therefore, we used measures closely 
related to each of Hill's first three diversity orders: zero 
(number of alleles, NA; allelic richness, _R), unity (Shan- 
non's Index, H), and two (Hardy-Weinberg expected 
heterozygosity, H E ). To calculate NA, R, and H E for each 
sample, we used FSTAT, version 2.9.3.2 (Goudet 1995, 
2002), and S H was calculated using GenAlEx. For greatest 
utility in future comparisons, we also convert diversity 
orders 1 and 2 into their effective number equivalents, 
which avoid many well-known problems of diversity mea- 
sures (Jost et al. 2010; Leinster and Cobbold 2011). The 
respective effective numbers equivalents are 'D^d,;,, = 
2A S H, and 2 D = 1/(1-H E ). 

We used nonparametric Mann-Whitney [/-tests to 
compare within-population diversity levels between the 
samples identified as sources for Australian introductions 
of both species because these data could not be made 
normal by transformation. Three approaches were used to 
assess diversity between populations. Pairwise F sx values 
were calculated in Arlequin for comparison with other 
studies that quote this measure. Pairwise values for Shan- 
non's mutual information index ( S H UA ) were calculated 
in GenAlEx. Compared with f sx , mutual information is 
known to be more robust to a wide range of population 
sizes and dispersal rates (Sherwin et al. 2006; Dewar et al. 
2011); additionally, the mutual information index can be 
converted to a numbers equivalent ('Db^ee,,), which 
avoids some serious problems that occur with other 
between-population measures (Jost et al. 2010). 

We surveyed the literature regarding H E measured from 
polymorphic microsatellite data in species from both fam- 
ilies containing our study taxa, Asteraceae and Caryo- 
phyllaceae, to determine whether H E estimates generated 
from native populations in this study were congruent 
with those from other members of the same family. This 
search was conducted in Google Scholar using the family 
name as a search term in conjunction with the terms 
"microsatellite" and "heterozygosity" in August 2012. 
Where data were given for multiple populations within a 
study, a mean value of H E was used. We avoided includ- 
ing estimates generated from introduced ranges, those of 
populations suspected of hybridization, and those of 
cultivated populations. Then, we surveyed the literature 
for examples of species showing evolutionary change in 
their introduced range, where genetic diversity had been 
estimated in both the native and introduced ranges. This 
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search was conducted in Google Scholar in October 2012 
using the terms "introduced" and "heterozygosity" in 
conjunction with either "rapid evolution" or "contempo- 
rary evolution." Additionally, we included studies refer- 
enced in a review of genetic variation across native and 
introduced ranges (Dlugosch and Parker 2008a) showing 
evidence of morphological change in the introduced envi- 
ronment. We calculated the ratio of diversity found in 
the introduced range to that found in the native range 
( He)> which gives an estimate of diversity retained after 
introduction, assuming no changes in diversity have 
occurred in the introduced range. For this calculation, we 
only used estimates of H E generated from microsatellite 
data because the absolute values of diversity estimates 
differ according to the marker used, and we wanted to 
directly compare these results to those generated in the 
current study. 

Results 

Evidence of morphological change 

Petrorhagia nanteuilii in the introduced range showed a 
significant increase in log 10 height through time that was 
not driven by region sampled (weighted general linear 
model including a term for region; R 2 = 0.06; 



iVeari,i76 = 5.06; P region = 0.13; P = 0.026; Fig. 3). 
Although the predictive power of this relationship is low, 
the magnitude of change is high, with average plant height 
increasing by almost 30% between 1882 and 1998. 

In the native range of A.populifolia, log 10 plant height 
was unchanged across the period of this study (weighted 
general linear model; R 2 = 0.003; -F ye ari,s2 = 0.15; P year = 
0.70; Fig. Sla). In the native range of P. nanteuilii, log 10 
plant height decreased through time (weighted general 
linear model; R 2 = 0.05; F yeltrlt g 5 = 3.98; P year = 0.049; 
Fig. Sib). 

Microsatellite markers 

We found no evidence for departures from Hardy-Wein- 
berg and linkage equilibrium in the microsatellite data for 
A. populifolia. Two of the twelve loci developed for P. nan- 
teuilii (Pna06 and Pnal6) significantly deviated from 
Hardy-Weinberg equilibrium and were excluded from 
downstream analyses. The remaining ten loci showed no sig- 
nificant departures from equilibrium. Petrorhagia nanteuilii 
has previously been reported to be tetraploid (Thomas and 
Murray 1981). Although we found no evidence of tetra- 
ploidy in the microsatellite data presented here, it is possible 
that in allotetraploid species, only a single parental genome 
may be amplified from any pair of primers. 
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Table 1. Estimation of diversity within native and introduced populations of Arctotheca populifolia and Petrorhagia nanteuilii including measures 
across three diversity orders: zero (allelic richness, ft), unity (Shannon index, S H; and the effective numbers equivalent, 1 D wit hi n ), and two (Hardy- 
Weinberg expected heterozygosity, H E ; and the effective numbers equivalent, 2 D). For both species, genetic differentiation within the native and 
introduced ranges was calculated using F ST and Shannon's mutual information index ( 5 H UA ), and the numerical equivalent of 5 H UA ( 1 Db e tween)- 



Statistic 



Arctotheca populifolia (A, 



Native 



Introduced 



Arctotheca populifolia (A 2 



Native 



Introduced 



Petrorhagia nanteuilii (B 2 



Native 



Introduced 



Mean ft (range) 

Mean 5 H (range) 

Mean 1 D withln (range) 

Mean H E (range) 

Mean 2 D (range) 

Mean Pairwise F ST (range) 

Mean Pairwise S H UA (range) 

Mean 1 D be tween (range) 



2.2 (2.0-2.6) 
0.79 (0.57-0.96) 
1.73 (1.48-1.94) 
0.34 (0.23-0.43) 
1.53 (1.31-1.74) 
0.33 (0.09-0.56) 
0.23 (0.08-0.41) 
1.18 (1.06-1.33) 



1.0 (1.0-1.1) 
0.08 (0.0-0.30) 
1.06 (1.0-1.23) 

<0.01 (0.0-0.01) 

1.01 (1.0-1.01) 
<0.01 (0.0-0.01) 
<0.01 (0.0-0.01) 

1.00 1 



1.3 (1.0-1.8) 
0.18 (0.0-0.46) 
1.14 (1.0-1.38) 
0.07 (0.0-0.18) 
1.09 (1.0-1.22) 
0.40 (0.19-0.59) 
0.09 (0.01-0.16) 
1.06 (1.01-1.12) 



1.1 (1.0-1.2) 
0.04 (0.0-0.09) 
1.03 (1.0-1.06) 
0.01 (0.0-0.03) 
1.01 (1.0-1.03) 
0.05 (0.02-0.08) 
0.01 (0.01-0.002) 
1.01 (1.00-1.01) 



0.09 (0.08-0.11) 
1.07 (1.05-1.08) 
0.04 (0.02-0.05) 
1.04 (1.02-1.06) 
0.56 (0.18-0.77) 
0.33 (0.27-0.41) 
1.25 (1.21-1.33) 



1.4 (1.3-1.4) 
0.33 (0.29-0.38) 
1.26 (1.22-1.30) 
0.13 (0.11-0.15) 
1.15 (1.12-1.18) 
0.52 2 
0.26 2 
1.20 2 



1 All measures equal mean. 

2 Denotes a single pairwise comparison. 



Population structure 

The Structure analysis of A. populifolia suggested the 
presence of two genetic groups (Fig. 4). One group 
included native samples extending from the western edge 
of the range in South Africa to Kenton on Sea as well as 
samples from south-eastern Australia (Fig. 4, group A^. 
The second group contained native samples from the 
eastern edge of the range to Kei Mouth and also included 
samples from Western Australia (Fig. 4, group A 2 ). Two 
genetic groups were identified in P. nanteuilii: One group 
consisted of most sampling localities in Spain and France 
(Fig. 4, group B^ and a second group (Fig. 4, group B 2 ) 
contained samples from two localities in Spain (Plasencia 
del Monte, PM; and Girona, GI), Corsica, Sardinia, and 

10 i 
9 - 
8 - 
7 - 




I I I 

0.1 0.2 0.3 0.4 



Australia. The samples from Corsica, Sardinia, and Aus- 
tralia had membership proportions for group B 2 in excess 
of 0.97, whereas those from Spain were lower (PM, B 2 
membership proportion = 0.61; GI, B 2 membership pro- 
portion = 0.80). 

Principal coordinate analysis plots were generally con- 
cordant with Structure results. Arctotheca populifolia 
samples from group A l (Fig. 5A, diamonds) formed two 
distinct clusters representing native and introduced popu- 
lations, respectively. A^introduced samples were more 
closely related to Aj native samples than any A 2 samples 
(Fig. 5A, circles). Within the A 2 group, the Kei Mouth 
sample was separated from all other samples. The remain- 
ing native A 2 samples were clustered with introduced A 2 
samples. Petrorhagia nanteuilii samples from group B! 

0 Asteraceae ■ Cctryophyllaceae 




Expected heterozygosity 

Figure 6. Native range estimates of expected heterozygosity from microsatellite data in Asteraceae (striped) and Caryophyllaceae (solid) families 
(see Table S1, for details). Arrows indicate the level of heterozygosity found in the native ranges of the species used in the present study 
(calculated from all sites sampled in the native range of each species). 
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Table 2. Species showing evidence of evolutionary change in introduced environments and for which genetic diversity was measured in native 
and introduced populations. The ratio of genetic diversity in the introduced range to the native range is given, and the direction of change is 
given (Trend). The statistics used to calculate diversity included allelic richness (/?), expected heterozygosity (H E ) and genet richness (GR). Effective 
number equivalents (E) have been calculated for native/introduced diversity. 

Genetic Diversity 

Species Native Introduced Ratio Trend Statistic E Reference 
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'These estimates are from primary introductions, which were the sources of secondary introductions ("introduced" values for these species). 
2 Mean includes all species having H E estimates. 



(Fig. 5B, diamonds) were clustered together. The two 
Spanish samples that had lower membership proportions 
for group B 2 in the Structure analysis were clustered with 
samples in group B x . Group B 2 included the introduced 
sample from Hampton, which was clustered with the 
native sample from Corsica, and the introduced sample 
from Albury, which clustered with the native sample from 
Sardinia (Fig. 5B, circles). 

Genetic diversity 

In the introduction to eastern Australia, native samples of 
A. populifolia had significantly higher H E (Table 1) than 
did introduced samples (A[ introduction: Mann-Whitney 
17, P < 0.01), while in the introduction to Western 
Australia, H E was not different between native and intro- 



duced samples (A 2 introduction: Mann-Whitney U, 
P = 0.70). Diversity was low within the second introduc- 
tion: We found a single genotype across all seven loci in 
one A. populifolia sample from the northeastern extreme 
of the native range (Umlalazi) and three samples from 
Australia (Beachport, Narooma and Wairo Beach). 
Assuming a single introduction of P. nanteuilii to Austra- 
lia, the estimated H E in the native samples was not differ- 
ent to that of introduced samples (B 2 introduction: 
Mann-Whitney U, P = 0.33). Values of R were similarly 
low in both species (Table 1). Despite the low values of 
genetic diversity, we detected that within populations, the 
total number of alleles detected for each species was not 
particularly low; 36 alleles were detected across the seven 
loci used for A. populifolia and 50 alleles across the ten 
loci used for downstream analysis in P. nanteuilii. This 
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highlights the strong genetic diversity found between sam- 
ples across the native ranges of both species (Table 1; F ST , 
0.33-0.56; S H VA , 0.09-0.33). 

We found microsatellite estimates of H E in 28 Asteraceae 
species from 28 different genera ranging from 0.22 to 0.88 
(mean H E = 0.58, Fig. 6 and, Table S2). Of the approxi- 
mately 2200 species in Caryophyllaceae (Schweingruber 
et al. 2011), microsatellite data exist for only eight species 
from five genera. We used expected heterozygosity mea- 
sures from all eight of these species ranging from 0.35 to 
0.93 (mean H E = 0.64, Fig. 6 and Table S2). In this study, 
values of H E across the native ranges of A. populifolia and 
P. nanteuilii were lower than almost all values we found for 
confamilials (Fig. 6; A. populifolia H E = 0.26; P. nanteuilii 
H E = 0.17) and both were significantly lower than the 
mean of these values for each family (one-sample Mests: A. 
populifolia, P < 0.001; P. nanteuilii, P = 0.001). 

We found 19 examples in the literature of species that 
had demonstrated evolutionary change in their intro- 
duced environment and where genetic diversity had been 
estimated in both the native and introduced range. Eight 
of these studies used allozymes, half of which found 
higher genetic diversity in the native range (Table 2). Of 
the eleven studies that used microsatellites, only one esti- 
mated a higher genetic diversity in the introduced range 
as compared to the native range. Ten of the microsatellite 
studies could be directly compared to our results (i.e., 
they provided estimates of H E ). In these studies, the ratio 
of diversity in the introduced range to the diversity in the 
native range was an average of 0.81 ( R H E range: 0.30-1.22, 
Table 2). In the present study, more diversity was found in 
the introduced range than the native source population of 
P. nanteuilii ( R H E = 3.25). However, R H E for A. populifolia 
was much lower than estimates found in the literature 
search (A t introduction: R H E = 0.01, A 2 introduction: 
R H E = 0.14). In fact, both introductions of A. populifolia 
had the lowest ratios of any example we found. 

Discussion 

Genetic diversity has been demonstrated to be positively 
correlated with invasion success (Crawford and Whitney 
2010), and standing genetic variation is believed to be 
important to invasive species' ability to adapt to novel 
environments (Barrett and Schluter 2008). However, it is 
becoming clear that introduced populations with very low 
neutral genetic diversity are sometimes successful invaders 
(Ren et al. 2005; Mergeay et al. 2006; Zimmermann et al. 
2010) and have the ability to adapt to their new environ- 
ments (Dlugosch and Parker 2008b; Harris et al. 2012). 
Here, we provide two examples of species that have estab- 
lished, spread, and adapted to the environment in their 
introduced range in Australia (Buswell et al. 2011 and 



Fig. 3), yet have significantly lower genetic diversity in 
their native ranges than do confamilials. Interestingly, for 
both of these species, the changes we identified in the 
introduced ranges were not found in the native ranges, 
ruling out the possibility that global processes are driving 
these changes. 

Despite the fact that the two separate introductions of 
A. populifolia described here had comparatively low levels 
of genetic diversity in the native source populations, these 
two introductions represent a larger percentage loss of 
genetic diversity than found in any introduction identified 
in our review of species, showing substantial morphologi- 
cal change in the introduced range. Similarly, introduced 
populations of invasive Japanese knotweed (Fallopia 
species complex) harbored very low genetic diversity at 
Amplified Fragment Length Polymorphism markers 
despite displaying significantly different phenotypes in a 
common garden setting (Richards et al. 2008). While 
greater levels of diversity may increase the likelihood of 
invasion success (Crawford and Whitney 2010; Jones and 
Gomulkiewicz 2012), it is clear that some introduced 
species, such as those discussed here, are able to become 
invasive and adapt to their new environments with very 
little neutral genetic diversity. This has important 
management implications because it demonstrates that 
even introductions from very small numbers of individuals 
have the potential to become invasive. 

The review we have conducted specifically examines 
loss of genetic diversity at introduction in species where 
some evidence of adaptive change has been documented 
in the introduced environment. It would be useful to 
compare the associated change in genetic diversity in this 
group with that of a group of species which has been 
introduced but has shown no evidence of adaptation to 
novel environments; we might anticipate that the latter 
would show more loss if diversity is important to adap- 
tive potential. Unfortunately, there is a bias in reporting 
which makes this difficult. However, we can compare the 
results of our review (19% loss of GD in the introduced 
range) with that of Dlugosch and Parker (2008a), who 
found 22.6% loss of diversity at introduction irrespective 
of evidence of adaptive change (Dlugosch and Parker 
2008a) (two-sample f-test, P = 0.72). This suggests that 
neutral genetic diversity is not important to adaptive 
potential in introduced species. 

The ability of populations with low current "neutral" 
diversity to evolve could be due to (1) retention of greater 
adaptive than neutral genetic variation due to either chance 
or balancing selection on adaptive variation (Reed and 
Frankham 2001), (2) contributions of mutations to selec- 
tion response, especially when selection lasts for more than 
twenty generations (Frankham 1980b, 1983; Hill 1982a,b), 
(3) loss of neutral genetic diversity after much of the 
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adaptation has occurred, or (4) some combination of these. 
Substantial adaptive genetic changes can still occur in 
populations subject to bottlenecks, and there is often large 
variation among replicates (Frankham 1980a). We are 
unable to distinguish between these hypotheses, but the 
likelihood of contributions from mutations that arose after 
introduction increases as the level of neutral genetic diver- 
sity in the introduced population decreases. 

Recent research suggests that epigenetic modifications 
(DNA methylation) may also play an important role in 
invasion success. Using the invasive Japanese knotwood 
populations discussed above, Richards et al. (2012) dem- 
onstrated that although genetic diversity was extremely 
low, significant epigenetic differentiation occurred between 
sites, suggesting a possible nongenetic mechanism for 
adaptation. Theoretical work on epigenetic selection mod- 
els indicates that increased phenotypic change can occur in 
populations with no genetic variation as a result of epige- 
netic changes (Geoghegan and Spencer 2012). In fact, Liebl 
et al. (2013) found a negative relationship between genetic 
and epigenetic diversity in introduced populations of spar- 
rows (Passer domesticus) and speculated that epigenetic 
variation may provide a mechanism for adaptation over 
the short time scales relevant to invasions. 

Source population identification 

Our results highlight the importance of determining 
source populations prior to assessing changes in genetic 
diversity between native and introduced ranges. Few stud- 
ies investigating this topic have done this, but levels of 
genetic diversity can be very different across a species' 
native range, as we found with A. populifolia. Identifying 
the source of an invasion assures that observed differences 
between introduced and native populations are not the 
result of diversity within the native range and prevents 
actual differences from becoming obscured (Dlugosch and 
Parker 2008a). Similarly, for studies attempting to identify 
contemporary evolution in introduced species, it is vital 
that the correct source population is identified in the 
native range. 

Our analyses identified two genetic groups of A. populifo- 
lia. One group consisted of western South African native 
samples and eastern Australian introduced samples. PCoA 
indicated that the latter were most similar to native samples 
from the south coast of South Africa ranging from Muizen- 
berg to Kenton on Sea. The second group contained sam- 
ples from the east coast of South Africa and introduced 
samples ranging from Margaret River in Western Australia 
to Beachport in South Australia. The two genetic groups 
identified in South Africa correspond perfectly to the Cape 
Seashore Vegetation (Group Al) and the Subtropical 
Seashore Vegetation (Group A2) described by Mucina and 



Rutherford (2006). A single genotype across seven loci was 
found in all 30 individuals sampled at Umlalazi (eastern 
South Africa), and this genotype was found in every indi- 
vidual sampled in Beachport and 39 of the 43 individuals 
sampled in Western Australia. These results indicate two 
separate introductions to Australia. This is supported by 
morphological data indicating that two forms of A. populi- 
folia exist in Australia, both of which are found in South 
Africa (Heyligers 2007). The Victorian coastline separates 
the two morphological groups (Heyligers 2007) as well as 
the genetic groups found in this study. 

Samples within native populations of both species were 
highly differentiated. This is possibly due to our inten- 
tional selection of species with restricted ranges so that 
we could sample comprehensively across the native 
ranges. Although restricted ranges and high levels of pop- 
ulation differentiation can be caused by limited dispersal, 
there is independent evidence that there may be a small 
amount of long-distance dispersal in A. populifolia. Heyli- 
gers (2007) argued that the distribution of A. populifolia 
morphs in Australia could be explained by dispersal of 
achenes via coastal currents and that historical records of 
first appearance showed an eastward progression of this 
introduction from Western Australia to South Australia. 
Given this evidence, one might expect to find a more 
cosmopolitan distribution of this species, but to our 
knowledge, A. populifolia is only found in southern Africa 
and Australia. 

The Structure analysis of P. nanteuilii indicated that 
two genetic groups exist in the native range, but only one 
of these was represented within Australia. Mediterranean 
island samples from Corsica and Sardinia were most simi- 
lar to samples from Australia, supporting a single source 
for this introduction. However, while the PCoA sup- 
ported the membership of group B 1 determined in Struc- 
ture (Fig. 5B), the samples contained within group B 2 
were not well clustered. In fact, the PCoA indicated that 
the introduced sample from Hampton was closely related 
to the native sample from Corsica, whereas the intro- 
duced sample from Albury was closely related to the 
native sample from Sardinia. This raises the possibility 
that two introductions of this species into Australia may 
have occurred. 

Comparisons of genetic diversity between 
introduced and native ranges 

Low genetic diversity in introduced populations can 
reflect a genetically impoverished source (Voss et al. 
2012). Our identification of the source populations for 
the introductions discussed here allows us to confirm that 
the sources were genetically impoverished. Within native 
and introduced populations of A. populifolia and P. 
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nanteuilii, values of H E were considerably lower than 
found in other species within these families (Table 1; 
Fig. 6; Table S2). One native sample of A. populifolia dis- 
played no gene diversity (i.e., had a single genotype across 
seven loci). Variation was even lower in the introduced 
range of A. populifolia than the already low variation in 
the native range, suggestive of a small number of foun- 
ders. 

Despite low levels of within-population neutral genetic 
diversity, both of these species are widespread in their 
introduced range in Australia. Similarly, Hardesty et al. 
(2012) found extremely low levels of neutral genetic vari- 
ation in highly successful introductions of Miconia calves- 
certs (mean H E = 0.07). However, because the ability to 
evolve in response to a novel environment may depend 
on the level of adaptive variation present, the relevance of 
neutral markers variation has been questioned (Reed and 
Frankham 2001). Although a positive correlation has been 
reported between quantitative trait (Q ST ) and microsatel- 
lite variation (Pst)> quantitative variation is usually higher 
(Merila and Crnokrak 2001; Leinonen et al. 2008) and 
the ability to predict Q ST increases with increasing values 
of F ST (Leinonen et al. 2008). Population bottlenecks are 
predicted to reduce additive genetic variance (Wright 
1951; Chakraborty and Nei 1982), but in some circum- 
stances, such populations may experience an increase in 
additive genetic variance for traits with at least some non- 
additive genetic variation (Willis and Orr 1993; Wang 
et al. 1998; Willi et al. 2006). While this shift in additive 
genetic variance may not always result in an increased 
ability to adapt to novel selection pressures (van Heerwa- 
arden et al. 2008), the combination of adaptation with 
low diversity at neutral markers indicating the presence of 
a bottleneck has been identified here and in other studies 
(Koskinen et al. 2002; Yonekura et al. 2007; Dlugosch 
and Parker 2008b). Frankham et al. (1999) showed that 
the effects of population bottlenecks on ability to evolve 
in response to environmental change closely followed 
neutral expectations. Further, a number of data sets indi- 
cate that genetic variation involved in adaptation to new 
environments is approximately additive (de Oliveira and 
Cordeiro 1980; Frankham et al. 1999), in contrast to the 
fitness variation in the environment to which populations 
have been adapted long-term, where there is usually a 
predominance of nonadditive variation (i.e., the occur- 
rence of increased additive genetic variation in bottlenec- 
ked populations probably does not apply to populations 
adapting to new environments). 

Finally, we speculate on whether intrinsic characteristics 
of the species might be affecting the diversity. Finding a 
genetically monomorphic sample in the native range may 
suggest alternate forms of reproduction across the species' 
distribution. When Roman and Darling (2007) examined 



successful introductions having decreased genetic diversity 
in the introduced range, 63% had reproductive abilities 
other than those involving sexual recombination. It is 
possible that A. populifolia may have the ability to repro- 
duce via apomixis, spread vegetatively, or self- fertilize, 
although none of these reproductive mechanisms have 
been reported in this species. Arctotheca populifolia has 
previously been reported as diploid (Norlindh 1967), and 
our investigations of ploidy in both native and introduced 
samples support this (see Appendix S2). This suggests 
that apomixis is an unlikely explanation for the observed 
genetic pattern because, among Asteraceae, this reproduc- 
tive mechanism is normally only found in polyploids 
(Noyes 2007). Vegetative reproduction has been reported 
in the congener A. calendula (Bossard 2000), but does not 
explain the biogeographic patterns of A. populifolia 
described here. Baker's Law (Stebbins 1957) states that 
self-fertilization should provide an advantage to coloniz- 
ing populations (Baker 1955) and could explain the 
genetic patterns we have identified in eastern South Africa 
and Australian introduced populations. Baker's Law is 
supported by three findings: i) increased frequency of 
self-compatible species on islands (Barrett et al. 1996), ii) 
species capable of autonomous seed production had larger 
invasive ranges (van Kleunen and Johnson 2007), and iii) 
species naturalized outside of their native range are more 
likely to self-fertilize than congeners only found in their 
native range (van Kleunen et al. 2008). 

Conclusion 

Considerable effort has been invested in identifying driv- 
ers of invasion success, including the importance of 
genetic diversity to invasiveness. While genetic diversity 
may be related to invasion success in some species (Craw- 
ford and Whitney 2010), increasingly, evidence suggests 
that genetic diversity is not essential to a species' ability 
to invade novel environments. Here, we have identified 
two species with low levels of neutral genetic variation in 
both their native and introduced ranges, which appear to 
have adapted and spread in their introduced range. 
Recent empirical evidence and simulations suggest that a 
number of factors influence the relationship between 
genetic diversity and invasion success and that complexi- 
ties such as competitive interactions and diversity of the 
native community are likely to be important (Chang and 
Smith 2012; Hovick et al. 2012; Jones and Gomulkiewicz 
2012). Further, it appears that epigenetic modifications 
may play a role in facilitating invasion success immedi- 
ately following invasion, although this idea has not yet 
been rigorously tested. In combination, these results sug- 
gest that genetic diversity measures alone are inadequate 
predictors of invasion success. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Appendix SI. Theory of adaptive genetic change in novel 
environments. 

Appendix S2. Chromosome investigations of Arctotheca 
populifolia. 

Table SI. Characterization of microsatellite loci in Arctot- 
heca populifolia (N = 348) and Petrorhagia nanteuilii 
(N = 377) including locus name, GenBank accession 
number, primer sequences, repeat motif, number of alleles 
and allele size range. 

Table S2. Expected heterozygosity (H E ) estimates from 
microsatellite data for species within the Asteraceae and 
Caryophyllaceae families. Effective number equivalents 



( 2 D) have been calculated. For Asteraceae, estimates were 
included for one species per genus {N = 28). For Caryo- 
phyllaceae, all estimates identified in the literature were 
included (N = 8). The number of loci (L), samples (S) 
and total number of individuals (I) are given for each 
study. 

Figure SI. (A,B) Plant height (logi 0 transformed) of Arc- 
totheca populifolia (A) and Petrorhagia nanteuilii (B) mea- 
sured from herbarium specimens sampled in the native 
range between 1891-2003 and 1848-1985, respectively. 
For A. populifolia, values do not change significantly over 
this time period (weighted general linear model; 
R 2 = 0.003; F y eari,52 = 0.15; P = 0.70). For P. nanteuilii, 
plant height decreased through time (weighted general 
linear model; R 2 = 0.05; F yeaih&5 = 3.98; P year = 0.049). 
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